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ABSTRACT 


Submitted to the Department of Naval Architecture and Marine Engineering 
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Master of Science Degree in Electrical Engineering and the Professional- 
Degree, Naval Engineer* 

The Deep Submergence Rescue Vehicle is being built vith an activated 
roll stabilization system. The effectiveness of this system is limited 
by pump saturation. This investigation is concerned with improving the 
roll stabilization of the DSRV through the use of a passive tank stabi¬ 
lizer. A mathematical model for the tank system is developed and adapted 
to the DSRV model. The passive system is extended to include two modes 
of operation, each representing a separate tank frequency, vrhich may be 
selected by the operation of a valve. The passive system and the dual 
mode system are simulated on an analog computer to deteronine their res¬ 
ponse to transient disturbances. 

The results show that a passive stabilizer of reasona.ble design could 
improve the roll stabilization of the vehicle. Furthermore, it is con¬ 
cluded that stabilization of transient disturbances can be improved by 
dual mode operation if a tank frequency at least three times greater than 
the natural frequency of the vehicle can be obtained for one of the two 
tank modes. 


Thesis Supervisor: George C. Newton 

Title; Professor of Electrical Engineering 



















NOMENCLATURE 


A) 

Variables and Forcing Rmctions 

x,y,z 

Coordinate directions, vehicle body reference frame: x forward, 

y to starboard, z down 

• U,V,W 

Fluid velocities, inertial reference 

4> 

Roll angle, positive to starboard 

G 

Pitch angle, positive bow up 


Yaw angle, positive bow to starboard 

s 

Curvilinear fluid displacement along centerline of tank 

K 

Roll-producing moment applied to vehicle 


Moment applied to tank fluid 

K 

V 

Moment applied to vehicle resulting from tank fluid motion 


B) Parameters Related to Tank Geometry 


A(S) 

Tank cross-sectional area as a function of S 

A 

Pipe cross-sectional area 

Aj, 

Side tank cross-sectional area 

^Ave 

Proportionality constant relating force to fluid velocity 

e 

Tank fluid angle 

• 

d 

Perpendicular distance between roll axis and any po nt on tank 

centerline 

Dp 

Energy dissipated in fluid flow 

D 

P 

Damping term: a coefficient which includes effects of flow loss 

/g/L 

Natural frequency of fluid motion in tank system 

Y 

Angle between r and d. 

I' 

X 

Vehicle moment of inertia about roll axis modified by removal 
of passive system fluid 












IV 



Tank moment of inertia about roll axis 

K, 

A parameter reflecting conservation of fluid angular momentum 

L 

Tank effective length 

Jl - 

Pipe length 


Tank curvature constant 

V 

Tank moment about pitch axis 

n 

Number of connecting pipes 

R 

Side tank lever arm 

r 

Radius to any point on tank centerline 


Fluid kinetic energy 


Fluid potential energy 


Tank volume 

'^T 

Total fluid velocity 

w' 

Vehicle weight modified by removal of passive system fluid weight 


C ) Other Notation Used 


g 

Acceleration due to gravity 

h 

Differential height between two fluid levels 

«f 

Head loss in feet 

t 

Time in seconds 
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CHAPTER I 

GENERAL DISCUSSION OF RO L L STABILIZATION 

Introduction 

Roll stabilization of ships at sea has been an area of engineering 
interest for many years. Historically, the first recorded attempt to 
build a functional roll-stabilization system occuired during the latter 
years of the 19th century. Since that time there has been a steady 
technical effort directed at improving roll-stabilization techniques 
and equipment. Roll-stabilization has been recognized as a problem not 
unlike conventional, servo mechanism problems.^ Accordingly, more recent 
technical advances within the discipline of automatic control have been 
used successfully to overcome many of the control problems which plagued 
early systems. 

Currently, the most widely used roll-stabilization systems are gen¬ 
erally of the fin type or the moving fluid type. However, moving solid 
weights and gyroscopes have also been used in the past as the major 
element of stabilization systems. The use of fin type stabilizers is 
somewhat restricted, in that their effectiveness decreases rapidly with 
ship speed. However, fin stabilization systems have been used success¬ 
fully for many years aboard ships which operate primarily at or near 
their designed speed. Moving fluid systems (unlike fins) provide reliable 
stabilization at all speeds. These systems are grouped into one of two 
categories. A passive system is one which depends upon gravity and 
inertial effects alone for the transfer of fluid mass. An activated 
system employs some type of pumping device to force the transfer of fluid 
mass within the system. 

In either the active or the passive system, the principle of operation 
is the same. Energy is imparted to the ship system as a result of wind 

















and wave action or by other transient forces and moments such as those 
produced by ship control surfaces. The added energy is partially trans¬ 
ferred to the stabilization system's fluid mass. The fluid system is a 
conservative system, except for the energy dissipated due to flow losses. 
The resulting fluid motion develops forces and moments which oppose the 
ship motion and thus produce the stabilizing effect. The effectiveness 
of an anti-roll system is therefore dependent on the system's ability to 
transfer energy to the fluid, the fluid-to-ship motion phase relationship, 
and the energy dissipation characteristics of the fluid flow. Energy 
transfer to the fluid and the flow losses are primarily a fimction of the 
tank geometry. The phase relationship between fluid and ship motion de¬ 
pends primarily on tank geometry for passive systems. For active systems, 
the pump and control system responses, as well as the tank geometry, deter^ 
mind the- phase relationship. 

In general, activated systems are capable of greater roll reduction 
than are passive systems under similar conditions. However, as in most 
physical systems, a penalty must be paid for higher performance. In 
activated roll-stabilization systems, the penalty appears in the form of 
the weight and power requirements of the system's pump. 

The DSVR and Roll Stabilization 

The Deep Submergence Rescue Vehicle (DSVR) is being designed with 
control capability in six degrees of freedom. The synthesis and simula- 
tion of this system has been reported in the literature. Of primary 
interest in this investigation is the proposed DSVR list-control system, 
which, under cruising conditions (or hovering with zero list), would 
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effectively be an activated roll-stabilization system. List control in 
the propssed DSVR system is to be effected by pumping mercury between a 
set of three tanks. Appendix I illustrates the system model. 

Pertinent to the subject of roll stabilization of the DSVR are some 
of the physiceil characteristics of the vehicle. Basically, the DSVR will 
nave a body-of-revolution hull form with few significant appendages. The 
transverse metacenter is very near the axis of symmetry of the body of 
revolution, and as a result, the natural damping of the vehicle is prac¬ 
tically negligible. The lack of natural damping places a heavy burden on 
the list control system, because the system must then compensate for any 
and all disturbances in the roll mode. 

The DSVR, like many other deep submersible vehicles, is power and weight- 
limited. for this reason there is a physical limitation as to the size 
and capacity of the pxamp used in the list control system. Therefore, in 
operation, the system is subject to ptimp saturation which severely degrades 
the vehicle stabilization in the roll mode. Furthermore, because of 
nearly continuous pump operation, the power requirements of this system 
become an important factor in overall vehicle endurance. 

This investigation will be concerned with the adaptation of a passive 

roll stabilization system which will complement the existing DSVR roll- 

control system. A tank system similar to a conventional U-tube stabilizer 

2 

will be considered initially. This system will be extended to employ a 
valve arrangement which may be controlled in such a manner eis to utilize 
the effects of gravity and inertia to enhance the transfer of fluid mass 
within the system. The purpose of the valve is to increase the rate of 
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mass transfer and thus counteract the effects of pump saturation as 
experienced hy the existing DSVR roll-control system. The investigation 
will use as a "basis the DSVR geometry and vehicle model as reported in 
reference four and outlined in Appendix I. The modified model which will 
be assumed will not alter the exterior characteristics of the existing 


model. 












CHAPl’ER II 


MODELING THE PASSIVE TANK SYSTEM 


Assumptions 

It is assumed that a second tank system is designed and added to 
the existing DSRV vehicle. The existing roll-control system is left 
intact except that a portion of the existing mercury will be used in the 
second system. It is also assumed for the purposes of this investiga¬ 
tion that the total vehicle mass is unaltered by the addition of the 
passive system. The tank system will be assumed similar in design to 
the U-tube type shown in Fig. I. 


The Mathematical Model 

The general equations describing the fluid motion are developed 
using a Lagrangian approach. Three basic expressions are required to 
define the total energy of the fluid mass. Thus, expi’essions for the 
kinetic energy, the potential energy, and the energy dissipated due to 
fluid motion are required. Throughout the following development, trans¬ 
lational velocity of the vehicle and fluid in the x-direction is neglected. 

To define the kinetic energy of the fluid, an expression for its vel¬ 
ocity is required. A fluid particle has the following velocities in the y- 
and z- directions (see Fig. l); 


V = y + \i;X^ cos \f/ + S sin (3- (f)+y) ~ ^ sin (3 -<?>) (2.1) 


1 , ^ 


w = z _ 0 Xq cos 6 + S ^ 


cos (3-<j)+Y) - ^r cos (3-<j)) (2.2) 


z 

Since kinetic energy depends on the square of total velocity (v^ = v + w ), 


an expression for v^ is required. 














tz 


Figure I 

Velocity Components of a Fluid Particle Within a U-Tube Tank 
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ri = z^- 2zXfiecos 6 + cos^e + 2z S cos (3-(}) + y) 


- 2z4)r cos (B-4>) - 25 XqS^j cos (B-(|>+y)cos 


( 0 ) 


• + 29X (|)r cos 0 cos (3-<|)) + S 

0 


S^(^i?/^(g))^cos^ (3- <1> + Y) 
- 2|^^T/A(S)js4>r cos (3-4>) cos (3-4>+Y) 


+ ip^r^ cos^(3 -4>) + 


S^[^T/A(g)] sin^(3-4> +y) + sin^(3- ((>) 


2y^X^ cos (if)) + 2yS 


sin (3- <i> +Y) 


- 2y^r sin (3-^) + 2^X^s|^T/A(s) cosif) sin(3 - <f> + y) 


- 2'f)X^<|)r cos if) sin (3-<l>) 


- 2s|r[^> 


/A(S) 


sin (3-<f>+Y) sin (3-<|>) 


(2.3) 


Making the approximation that fluid velocity is constant at any cross- 
section within the tank system, the kinetic energy of such a cross-section 
may be expressed as follows; 


dT. 


= |^T/2j 


dm 


where dm = pA(S)ds 

Therefore, total fluid kinetic energy may be expressed as; 


= I 


2.k 
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Substituting Equation 2.3 and noting that the tank volume (v^) is j A(s)ds 
produces the following expression: 


~ 2^ {z^- 2zXgecos 9 + cos^ 9 + co3^ ip* 2yi>X^ cosHv^ 

- p4»(y + 'i'X^ cos l|;) I A(S)r sin(8-({»)ds 

- p4)(z - 6 X^ cos 9) |A(S)r cos(3-({>)ds 

+ pS A^(y + ifx^ cos ij>) |sin(3-(ti +y )ds 

+ pS A^(z-9Xq cos 9) |cos(3 - <(>+Y)<is 

+ —^— I pT /A(S)]‘^® “i" I + PAjiS^*! r cos Yds 


This expression may be simplified by using the approximations sin (j) = (|) 
and cos <{) = 1. These approximations are valid voider conditions of small 
roll angles which are to be expected in a stabilized vehicle. The first 
integral in Equation 2.5 thus becomes: 


A(s)r sin(3 - <|))ds 


vl 


A(S)r sin 3 ds + (p ( A(s)r cos 3 ds 


The first integral on the right above is recognized as the tank moment 
(M ) about the pitch axis. The second integral is, by inspection, zero 
because of the tank symmetry. Therefore; 


A(s)r sin 3ds = M, 


Ty 


2.6 


Similarly: 


j A(S)r cos(3 -4*)ds = | A(S)r cos 3 ds - <|) |A(s)r sin 3 ds = 


2.7 


The third and fourth integrals in expression 2.5 are seen (see Fig, l) 


to reduce as follows: 













I cosO - <}) “ 


sin(3 - <j> +>0*^3 “ 


I cos(3 + Y)ds - <j>| sin($+Y)ds 


0 - <}>2R 


I sin(g + Y)ds + <J)| cos(3 + y)^s 


2R + 0 


2.8a 


2.8b 


where R is the tank radius. 

Investigation of the three remaining integrals in expression 2.5 
reveals that these represent constants that are unique to the chosen geo¬ 
metry of the tank system. The following definitions are made: 


di 

["S^/ACS)] = 

L 

(Effective length) 

2.9 

1 

r cos Yds = 


(Curvature constant) 

2.10 

1 

pr^A(S) ds ^ 


(Tank moment of inertia) 

2.11 


Substituting equations 2.6 through 2.11 into Equation 2.5 yields the 
kinetic energy expression: 


Tf = (z^- 2 zXq0cos e + 0^Xg cos^6 + y^ + cos^iIj + 2y^X^ cos 

+ (y + 4)Xj^ cos ijj)(+ pSA^2R - p^M^) 

- (z - 6Xq cos 0)(4)2RpSA^ + p|<J)M^) 

+ pA^S^L + Jj pAj,m\ 2.12 

The expression for potential energy may be considered to have three components 
each of which may be calcxilated separately. With the tanks blocked (S = O), 
the potential er ergy of the fluid due to roll is given by the expression: 
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Letting cos (|) 


1 

The potential 


Letting sin 6 


The potential 


= 1 —^ this equation reduces to: 
pg I A(S)r sin g ds = 

energy contribution due to pitch (tanks blocked) is: 

= |peVe sin^e. 

2 

= 0 this becomes: 

Vf = |pEV^X00^ 2. lit 

2 

energy contribution due to fluid motion in the tanks is given by: 


= pg^j S(S cos (J)- 2R sin<|)) 

^ 3 

As previously stated, the total potential energy is the Siam of the three com¬ 
ponents. 


V, = V ♦ Vf - V 

12 3 


^ Ajy + I PgV^Xge^ + PgA^ S(S - 2R<(>) 


2.15 


The dissipative energy expression will be postulated to be a fimction of the 
cube of the fluid velocity relative to the tank system. 



The function is not yet specified. 


2.16 














The generalized coordinates vhich are of principal interest in this 

$ 

development are (p and S . The two Lagrange equations of interest are, 
therefore: 


d 1 


1 '■'f. 




dt 1 

[srj 

■ 34. 

d(j) 

^f 

2,17 


Kl 

3T 

3Vf 

3D 


dt 

F* 

"33 

3S 

=0 

2.18 


It should be noted that Equation 2.17 is a moment equation and that 
represents the moment (in the <|) direction) applied to the fluid of the 
tank system. Equation 2.l8 is a force equation representing the force 
applied to the fluid in the S direction. 

Using T^, V^, and as defined by Equations 2.12, 1.15 and 2.l6, 
respectively, in Equations 2.17 and 2.l8 produces the following two basic 
system equations: 

+ PgHj,y4' - P\M^S - 2pgRA^S + (z - SXq cos 6)2RpA^S 
- P^'’5.y(y + cos > + 4)z - <t.eXg cos e) = 2.19 

2pA^,LS + + 2pgA^S - 2pgRA^<f) - pApM^cf) 

•• 

+ 2pA^E(y + ij;X^ cos ij; - z$ + |eXg cose - <j)Z + <l)eXQ cos 6)= 0 

2.20 

These equations will be simplified by considering the terms containing 
the variables z, z, y, 0, 6 and ip to be disturbance inputs. The value 
of any one of these variables would be small during any reasonable vehicle 
maneuver and a first-order approximation would be to assume that all terms 
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containing the product of two vziriables are negligibly small. However, 

• 2 

since the S tem in Equation 2.20 will be retained as the tank damp-- 
ing term, it would be inconsistent to set all product terms to zero. 
Alternatively, it will simply be assumed that the disturbance inputs may 
be individually set to zero so that attention may be focused on the rela¬ 
tionship between vehicle roll and fluid displacement. In addition to the 
simplifying assumptions. Equation 2.19 is altered by noting that any 
moment applied to the tank fluid must be opposed by an equal moment acting 
on the tank. Thus, Equations 2.19 aJ^d 2.20 become, respectively: 

J5.4)‘+ pA^S - 2pgRA^S = - \ 2.21 

*i •• 

2 pA^LS + S^+ 2pgA^S - A^^<t> - 2pgRA^<}> = 0 2.22 

The forcing term in Equation 2.2l(K^) is the moment applied to the vehicle 
as a result of fluid motion within the tank system. 

Adapting the Tank Model to the DSRV 

The simplified equation of motion in roll for the basic DSRV (see 
Appendix I) is; 

• • 

(l - K)4> = -wz<t> + K +K |ul^ 2.23 

X p ^ g S plW C-.C.J 

The coefficients of Equation 2.23 must be modified to reflect the change 
in the amount of mercury assigned to the basic vehicle model. The effec¬ 
ted terras are I and w which will be designated I' and w* to denote the 

X X 

modification. The total moment (Kg) applied to the vehicle is considered 
to be the sum of two components 

Kg = K^ + K 2,2k 

Combining equations 2.22 and 2.23 by using 2.2U yields the basic moment 

equation for the vehicle modified by the addition of the passive tank 









system. 




.(I* 

X 


- - Kpi^j |ul^+{V7^ + pgM^}(j) 


- pA^'^S - 2pgRA^,S = K 


2.25 


Equation 2.25 together with Equation2.22 defines the relationship he- 
tween ship roll and fluid displacement. 

Recall that Eq. 2.22 is a force equation where each term represents 
a force acting on the fluid. The term represents the friction force 

due to fluid flow. Since this force always opposes fluid motion, it can 
he approximated by the following expression; 

Kjjp |S|S 2.26 

The approximate proportionality constant is evaluated in Appendix II. 
Substituting Eq. 2.26 in Eq. 2.22 and rearranging yields: 

S = 


I* 2 pgA^ 
Similarly, Equation 2.25 becomes: 


ls!s - s + I; 5 + R({>| 


2g 


2.27 


4> = 


or 


I' - K* + 

< X p T/ 


+ Kp|^||u|^ - (w':^+ pgM^)({) + pAjM^S 
2 pgRA^S + kJ 


where: 


4) = + A^-A(|) + AS + AS + A K 

X 2 3 4 5 


■^s ^d' - K,+ J ) 

x p T 


A = K , Jul A 
1 p |ujl I 5 


A = (w'z 




A = 




2.28 

2.29 

2.33 

2.31 

2.32 



iHiniiiriiii 















lit 


A = 2 pgRA, A 

.4 X 3 


2.33 


Equations 2.27 and 2.28 are in a form suitable for simulation on an 
analog computer. Furthemore, in this form the significance of the 
various tank parameters is readily seen. The coefficient g/L in Eq. 2.27 
will be defined as the square of the tank natural frequency (as indeed 
it would be, if the equation were linear). Note also that, from Equa¬ 


tions 2.25 and 2.31, the coefficient A is the squared natural fre- 

2 


quency of the vehicle modified by the tank system. The coefficient 
Kjjp/2 gA^ will be referred to as the tank damping term. Tank frequency 
and damping are discussed in the next chapter in connection with the 
selection of the tank geometry. 
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CHAPTER III 


SELECTION OF TANK GEOMETRY MD PARAMETERS 

As previously stated, the passive tank system design is to he com¬ 
patible vith existing DSRV external geometry. Some of the basic design 
considerations are fairly self-evident. For exaraple, to result in maxi¬ 
mum effectiveness, the system should utilize the maximum ship beam, and 
locating the tanks at or near the ship’s yaw axis will minimize the inter¬ 
coupling effects between fluid motion and yaw. It is assiamed that the 
tank can be loce,ted on a transverse plane that includes the yaw axis. 

It is also assiuned that no portion of the system will be more than 3.5 
feet from the vehicle's longitudinal eu;is. Tlie latter constraint is 
required to keep the tank system within the exterior skin of the vehicle. 

A basic parameter in tank selection is the volume of fluid the sys- 
ten is to use. Previous experience has indicated two fundamental design 
points for tank selection.^ The first such point is that a total fluid 
weight of between 1/2^ and 1-1/4^ of total ship displacement is suffici¬ 
ent to produce effective roll stabilization if the tanks are properly 
arranged. Design point two is tha,t for a steady list of 1®, the tank 
static moment should be at least 12^ of the moment to heel 1®. 

To satisfy the first point, it is assumed that 1.66 cutric feet of 
mercury (equal in weight to 1% of total vehicle displacement) will be 
used in the passive system. This, by assumption, reduces the amount of 
mercury in the existing D3VR roll-control system by 50^. Allowing for 100/S 
fluid transfer under extreme conditions, each side tank volume must then 
be at least 1.66 cu. ft. A tank area of 0.6 square feet was selected, on 
the basis of being a compromise between reduced tank lever arm and the 
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overall tank length. The tank is assumed to be round in cross-section, 
with a curved vertical axis so as to conform to the geometry of the DSVR. 
Figure II illustrates the geometry. The small tank area and the assumed 
free surface have a negligible effect on vehicle static stability. This 
assumed tank arrangement results in a static tank moment for 1® of heel 
equal to about 25^ of the moment to heel 1®, and is thus well within the 
guidelines of design point two. 

The selection of the geometry for the cross-connecting portion of the 
tank system becomes somewhat more involved. Since this portion of the 
system constitutes a major portion of the tank’s overall length, the shape 
and size of the cross-connect directly affects the tank parameters 
L, and J^. In addition, the tank damping term is primarily a function 

of the geometry of this section. This latter point results from the fact 
that the damping term is dependent upon fluid friction which, in turn, is 
dependent on surface area exposed to the fluid flow. 

' In selection of the cross-connect geometry, four factors were con¬ 
sidered. For the side tank arrangement shown in Fig. II, a maximum (case l) 

J 

and a minimum (case 2) overall pipe length were selected for investigation. 
Six pipe cross-sectional areas were considered, which correspond to inside 
diameters of 0.08, 0.09, 0.10, 0.11, 0.12 and 0.13 feet. These values were 
selected because, in the author's judgment, they represented realistic values 
for such an application. Another factor considered was the total number of 
connecting pipes. Calculations were made for the use of one through eight 
connecting pipes; however, pipe diameters were not mixed for any single cal¬ 
culation. The last factor considered in the tank parameter calculations was 
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Figure II 


Basic Tank Geometry 


































pipe roughness. The "basis for all damping term calculations was a Moody 
diagram (see reference page 15 ) which is a plot of friction factor (f) 
versus Reynolds number. The curves given for smooth pipe flow and fully 
rough flow were used in these cailculations to illustrate the range of the 
expected damping term. 

The first step in the calculations was to solve equations 2.6, 2.9 - 
2'. 11 for and respectively. This was done for both case 1 

and case 2 with n (the number of pipes) and (the pipe cross-sectional 
area) as variables. These calculations and results are docvmiented in 
Appendix II. 

The calculation of the damping term is based on a theoretical maximum 

fluid velocity for static conditions of maximum heel. For transient fluid 

flow from a tank xmder the conditions of no friction, no body forces and 

tank diameter much larger than pipe diameter, the velocity in the pipe is 

7 

given by the expression: 

V= / 2gh tanh (t/^2gh/2)l) 3.1 

The maximum velocity for the assumed geometry iinder conditions of 30® of 

.* 

heel is therefore l6 ft/sec. imder ideail flow conditions. This figure is 
used as a maximiim vailue in calculating the damping term. It should be 
noted that from Eq. 3-1 about four seconds, or about 35^ of the DSVR’s 
natural period of roll, is required to reach the maximum velocity. 

The force exerted on the fluid flow due to friction was calculated 
by using the equation: 

Pounds force = pgA^H^ = pA^ {1.5 + f(l/d)} V^/2 3.2 

In the above expression, 1/d is the pipe length to diameter ratio, f is 
the friction factor obtained from the Moody chart, and the 1.5 represents 













an approximation of the entrance and exit losses of the flow* The results 
of Equation 3.2 (for velocities between 1.0 and l6.0 feet/sec) were 
averaged to determine an approximate proportionality constant. The re¬ 


sult is the following: F = * Note that from Eq, 2.27 term 



Since the tank velocity and the pipe velocity are related by the ratio of 
the respective areas and the nviiuber of pipes ^ the force term may he written 
as: 



2 


3.3 


V 


From the above it is clear that: 



3.h 


The relationship between the square of tank frequency and the tank damp¬ 
ing term are shown for case 1 and case 2 in Figures III and IV respec¬ 
tively. 

The two geometry cases-considered were selected because it was felt 
they would establish realistic end points. Case ;two, because of the mini¬ 
mum effective tank length (single pipe), produced the minimvim tank period. 
Therefore, maximum mass-transfer rates are to be expected from tlis geometry. 
No attempt was made to optimize tank geometry. 

In a practical system, tank tuning could be accomplished by varying 
the amount of fluid in the system, by selective switching of various cross- 
connect length/number of pipes, or by using valves to throttle the flow in 
the cross-connect pipes. Yet another possibility for tank tuning is to 
have a secondary fluid on top of the tank fluid (for a mercury system, oil 
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coxild be used as the secondary fluid.) Timing could then be accom¬ 
plished by throttling the return flow of the secondary fluid. How¬ 
ever, any of these possibilities could be used only to decrease the tank 
frequency from the maximum possible from the case II geometry. 

The calculations for pipe damping should be recignized as approxi¬ 
mate rather than absolute values. The damping term results represent 
the range of values which could be expected from a physical system. 
However, in an actual system design, more accurate results could be 
obtained from physical model tests. 
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CHAPTER IV 
SYSTEM SlfmATIOH 

The Passive System 

* 

The first simulation attempts of this system were conducted under 
conditions which were calciilated to show the effectiveness of the 
conventional passive system. Equations 2.27 and 2 . 28 , which represent 
the passive system, were used as the basis of this simulation. 

The system coefficients as defined by Equations 2.29 through 2.33 
were evaluated on the basis of the previously selected tank geometry. 

The coefficient Aj has a very weak dependence on the tank moment of 
inertia, since is less than 1^ of the vehicle's moment of inertia. 
Therefore this coefficient was taken as a constant (1.62 x 10 ®Slug"*ft 
throughout the simulation (see Appendix I for the numerical evaluation). 
The coefficient as defined by Eq. 2.30 represents a damping term which 
is dependent on vehicle motion in the horizontal plane. A worst case 
was assumed where the vehicle would be in the hover mode ( u = 0) and 
therefore A^ would go to zero. The coefficient A^ as defined by Eq. 

2•31 was evaluated and found to be rather insensitive to variations of 
the tank parameter As noted in Chapter two, this coefficient rep¬ 

resents the natural frequency (squared) of the vehicle modified by the 
tank system. The evaluation indicates that within the two assumed geo¬ 
metry cases, the natural frequency of the vehicle is approximately a cons¬ 
tant. The values used for A are: 


* The analog computer simulation of this system was carried out on a GPS 
model 200T analog computer. The simulated systemfe output was recorded 
on a dual channel, model 60, Sanborn Strip Recorder. Both pieces of 
equipment are part of the facilities of the Electronic Systems Laboratory. 
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Case i: A^' = 

2 

. 32 4 Sec“^ 

4.1 


-2 


Case II: A = 

2 

.319 Sec 



as defined by Eq. 2.32, has a strong dependence on the tank curva¬ 
ture constant. The values used for this coefficient are: 

Case I; A = .OO 89 Feet 

3 

Case II: A, = .OOUl Feet"^ 

3 

The coefficient A which is defined by Eq. 2.33, is dependent upon the 

4 

side tank geometry. Since this was assumed constant for both Case I 
and Case II, the coefficient is a constant. The value used is: 

A = 0.0507 (Feet/Sec^ )”^ U.3 

4 

The coefficients of Eq. 2.27 are all dependent on tank geometry. The 
tank frequency (squared) and the tank daiaping term are tabulated in 
Tables II and IV respectively. The other two coefficients, M^/2g and 
R, are self-explanatory. 

Initial simulation of the single mode passive system was conducted 
to study the effectiveness of the passive tank system. It was assumed 
that the vehicle roll-control system was inoperative, so that the full 
effect of the passive system could be observed. The basic system computer- 
simiilation diagram is shown in Fig. V. Scale factors have not been inclu¬ 
ded for simplicityInitial conditions simulating transient disturbances 
of heel (({)) and rolling ((f)) were applied to the vehicle model and the 
resulting damping effect recorded. The system model was also-subjected to 
a series of ramp and sine wave inputs, simulating a variety of disturbance 
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Figure V 


Passive System Simulation Diagram 






































































moments applied to the vehicle. Sample results of this simulation aie 
shown in Figures VI through IX. The effect of pipe roughness on tank 
effectiveness was observed and a sample result is shown in Figure X. 

A comparison in effectiveness between the existing DSVR roll-control 
system and the passive tank stabilizer was made. The roll-loop model 
reported in reference four was used to simulate the existing system. The 
responses of this model and the passive tank model are compared in Figure 
XI for similar disturbance inputs. The variable W in Figure XI is the 
pumping rate of the list-control piunp. This figure illustrates the 
major deficiency (piimp saturation) of the existing system. 

The Dual Mode System 

The simulation of the dual mode system presented some additional 
problems in that the effect of the valve closure on both the fluid and 
the vehicle motion must be accounted for. The valve(s) in the cross- 

connect pipe{s) which is suddenly closed while the fluid is in motion 

• •• 

would have the effect of driving S and S to zero, while S would remain 
at the value reached at the time of closure. This assumes an ideal valve 
which is either fully open or fully closed. In addition, the momentum 

j 

of the moving fluid column before the valve is closed must be at least 
partially conserved when the valve is closed. This effect was accounted 
for by assuming that the angular momentum of the. vehicle after the V 2 dve 
closes is equal to the vehicle's angular momentum before the valve closes 
plus a portion of the angular momentum of the fluid column. This rela¬ 
tionship can bo shown by the following expression: 

K" ^b" Vi® 


where ^ , k , k are the vehicle's angular velocity after the valve 

& D g 1 

closes, the angular velocity before the valve closes, a geometry constant 
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FIGURE m: 

RESPONSES TO INITMKL VELOCITIES 
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FIGURE VII 

Case I and case II responses to initial heal angle 
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Figure VIII 

Case II Responses to Ramp Moments. 
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FIGURE IX 

Response to periodic disturbances 
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FIGURE XI 


eomparison between existing system and the passive system, 
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relating S to angular momentum, and a variable (0 j< jS which will 
reflect the amount of fluid momentum transferred to the vehicle. A simu¬ 
lation diagram for the valve action is shown in Figure XII. 

Operation of the circuit shown in Figure XII depends upon two basic 

0 

characteristics of the 200T computer. When no mode control is applied 
to an integrator, it will act as a svamning junction for the signals applied 
to the initial condition inputs. Thus, when the valve is open, the out¬ 
put of 1-5 is - (j) + k k S. When the valve is closed, 1-5 becomes an 
integrator and the output includes the effect of the momentum impulse 
k k S. In the S- channel, when the valve is closed the switches change 

O 

position and mode control is removed from I-U. Because the disconnected 
output terminals of the 200T switch units are grounded, this drives fe and 
S to zero; It also resets I-U so that when the valve is reopened, S starts 
from zero. By grounding the input to 1-3 when the valve is closed, this 
integrator goes into a hold mode and the S value present when the closure 
occurs is retained until the valve reopens. 

The dual mode condition was simulated by making step changes in the 

j 

parameters g/L, and k^. The system was assumed initially to be in mode 

I, which has a short period (large g/L). This would facilitate a large 
mass-transfer rate when the vehicle was initially subjected to a large 
disturbing moment. At a later time, the tank system is switched to mode 

II, which is tuned to the natural frequency of the vehicle. The effect 
of switching modes is similar to complete closure of the cross-connect 
pipes, but not as severe, because the fluid column would not be completely 
blocked. In a physical system, mode switching could be accomplished by 
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Figure XII 



Valve,Simulation Diagram 
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closure of some of the pipes in a multiple-pipe cross-connect arrange¬ 
ment, or simply hy throttling the valve(s) in the system. Figure XIII 
shows a typical mode-control switching circuit. Note that A-6 and A-7 
are required because the electronic switch must be driven by a low impedance 
source. 

The selection of the valve timing and the mode switching criteria 
was approached on a trial and. error basis. Eventually, a very simple 
switching law was chosen which proved to give satisfactory results. The 
law is illustrated in Figure XIV. The system responses for sample initial 
disturbances and mode tunings are shown in Figures XV, XVI and XVII. 

Figure XVIII illustrates the effect of k^. 


j 
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CHAPTER V 

DISCUSSION OF RESULTS MD CONCLUSIONS 

The Passive System 

The results of the simulation of the passive system make two very 
clear points. The effect of tank tuning on roll-stabilization is most 
obvious. In those cases where the vehicle is subjected to transient dis¬ 
turbances, the most effective stabilization occurs when the tank is tuned 
to the natural frequency of the vehicle. For a periodic forcing function, 
the tank must be tuned to the frequency of the disturbance. 

The second point revealed by the basic simulation is that fluid trans¬ 
fer rates of 100 - 200 pounds per second are possible in the passive system. 
(Recall that the maximum pumping rate for the existing DSRV list-control 
system is only 56 pounds per second.) Observe, for example, in Fig. VI, 
for the initial velocity disturbance of 0.18 radians per second, the max¬ 
imum fluid excursion is about 1.32 feet over a time interval of about seven 
seconds. On the average, this is 96 pounds per second. Note that this 
example is for a tank tuned slightly below the vehicle natural frequency. 

t 

When the tank is detuned to a higher frequency, the flow rate is increased. 

The comparisons shown in Fig. XI indicate the increased effectiveness 
gained by using the passive system in lieu of the existing activated system. 
Fig. XI indicates that pump saturation occurs even with relatively small 
disturbances. However, this sample result also indicates that the effec¬ 
tiveness of the passive system is relatively independent of the magnitude 
of the disturbance. 
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An interesting point is illustrated in Fig. X. Otserve that, for 
the same conditions of tuning the smooth pipe system is more effective 
in stabilizing the disturbance than is the rough pipe system. This would 
indicate that there is an optimum value of flow losses which would result 
in maximum damping effect. This conclusion is reached because, if the sys¬ 
tem were lossless, there would be no energy dissipated, and therefore no 
damping effect. If the pipe were infinitely rough, there would be no 
fluid motion and again no damping effect. Therefore there must be an 
optimum value between these two end points. 

The Dual Mode System 

The effectiveness of the dual mode system is illustrated in Figs. 

XV through XVIII. In general, the amount of increased roll reduction 
from the dual mode operation versus the simple passive system is not large. 
One could question whether the increased system complexity required for a 
dual mode system could be justified in terms of the net gains in overall 
stabilization. Furthermore, the results shown here are the responses to 
some rather idealized, transient disturbances. 

» 

Nevertheless, there eire some significant results of this study of 
a dual mode system. The responses shown in Figs. XV, XVI and XVII indi¬ 
cate that increasing the mode I frequency increases the system's damping. 
These results also indicate that the mode I frequency should be at least 
three times larger than the vehicle natural frequency, if significant gains 
in effectiveness are to be achieved. 

Further study of the results of the dual mode operation reveals that 
the valve action increases the phase difference between the fluid motion 
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and the vehicle motion. For example, in Fig. XVI, observe that the 
time between the peak positive value of (j) and the maximum negative 
excursion of S is about nine seconds for the dusil mode system, and less 
than seven seconds for the same sequence in either mode I or mode II 
alone. Increased phase difference is the principal, result achieved from 
an activated anti-roll system. The point is that the power required to 
operate a valve would be considerably less than that required by a p\unp. 

Conclusions 

The major conclusion to be drawn from this investigation is that the 
amount of roll stabilization of the DSRV could be significantly increased 
through the use of a passive tank system. Under conditions of moderate 
disturbances, fluid transfer rates well in excess of 100 poiinds per second 
are possible for reasonable tank configurations. It has been shown that a 
tuned passive system is more effective in stabilizing moderate transient 
disturbances than is the existing system which is subject to pump satura¬ 
tion. In addition, it has been shown that there is a slight increase in 
system transient damping when a dual mode passive system is used instead 
of the simple passive system. The benefit of the dual mode system appar¬ 
ently results from the increased phase difference between the fluid and the 
vehicle motion. 

Suggestions for Further Investigation 

It is recognized that the existing DSRV roll-control system was de¬ 
signed to meet attitude control criteria as well as to stabilize the vehicle 
in roll. Furthermore, it is realized that any proposals for increasing the 
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roll stabilization of the DSRV must account for the capability to hover 
at any angle of heel less than 1+5® in magnitude. A possible candidate 
system vhich employs a passive tank stabilizer and still meets the heel 
angle requirements would be a passive tank system mounted on a subassem¬ 
bly capable of being rotated about the vehicle's longitudinal axis. The 
current list-control system could be retained for attitude generation 
and fine control, or the rotating subassembly coxxld be weighted in such 
a manner as to produce the required moment when rotated. 

In reference to the use of a dual mode system, further study is 
required to determine if the performance of a passive system subjected 
to the disturbances of an irregular sea could be improved throxigh dual 
mode operation. 











Figure XIII 

Mode Switching Arrangement 

. 



Figure XIV 

Valve Timing and Mode Control 
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FIGURE XV 

Comparison between dual and single mode operation 
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FIGURE XVI 


Comparison bet-ween dual and single mode operation 
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FIGURE XVII 

Comparison between dual and single mode operation 
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FIGURE XVIII 


Dua3. node response showing effect of k 






































































































































































































































APPENDIX I 


THE DSRV MODEL AND PARAMETERS* 

Nine generalized coordinates and generalized velocities are suf¬ 
ficient to model the DSVR system. For the vehicle without a moving 
mass center (no mercury) the velocities u, v, w, p, q and r describe 
the vehicle motion. The six needed Lagrange equations are: 


d 

ST 

V 

dt 

3u 
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ST 

V 

dt 

3v 
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3T 
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dt 

9w 
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1-5 


1-6 


The material in this Appendix has been drawn from reference li, and is 
included here for information only. 















T is the total vehicle kinetic energy and x, y, z, K, M, and W are the 

V • 

forces and moments acting on the vehicle. 

Three additional equations are needed to describe the effect of mer¬ 
cury motion. It was assumed that these equations could be reduced to 
quasi-steady equations defining the center of gravity location from the 
integrals of the various pvimping rates. 

The vehicle kinetic energy is given by; 


T 

V 


My " 

2 G 
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,y are the vehicle mass, the center of gravity velo- 
^G 

city, the angular velocity vector about the body axes, and the inertial 


where M, V_, w and 

vi 


tensor about the center of gravity. By defining axes x’, y’, z’ parallel 
to the body axis and through the center of gravity, 



^x'x' ^x’y' ^x'z' 

^y*x' ^y'y* ^y*z' 

'x' ^z 'y' ^z ’ z' 


The components of the velocity vector are: 


1-8 


\ = [u + qZQ-ry^j + i^,^ v+rx^-pzg+y^, w + py^- qxg+ Zq] i _9 

The inertial tensor about the center of gravity may be related to the 
inertial tensor ^about the center of buoyancy by the parallel axis theorem. 
For example 


I = 

XX 


x'x’ 


+ M(y^ + z^) 


I-IO 


By modifying Equation 1-8 through the use of a set of equations similar to 
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I-IO and using 1-9» the expression for T^(Eq. 1-7) may be solved. 
Substituting this result into Equations I-l through 1-6 yields the full 
set of dynamical equations for the DSRV. From Equation 1-4, the results 
for the roll mode are: 

K = I p + I (q-pr) + I (r+pq) + I (q^- r^) + (l - I )qr 
xx^ xy ^ ^ xz yz zz yy 

- MZg(v -pw +ru +rxQ ty^) + My^(w +pv -qu -qXg+ z^) 

+I p+I +I r+ Mx_(z„r + y_q) I-ll 

xx^ xy xz G' G ''G^ 

Equation I-ll was simplified by making the following assumptions: 

a) Because of xz plane symmetry, I =I =I =I =0 

^ ’ xy yx yz zy 

b) Assume almost fore-aft symmetry in weight distribution 

1 = 1=0 
xz zx 

c) Heglect velocity arid acceleration terms of the center of 

gravity: x^ = y^ = z^ = x^ = y^ = z^ = 0 

Note that because of this assiimption, I = i =1 =0. 

^ ’ XX yy zz 

t 

As a result of the simplifying assumption. Equation I-ll becomes: 

^xx^ ^^zz“^yy^*^^ " Mz^(v-pw+ru) + Myg(w+pw-pu) = K 1-12 

The forcing term in the above equation has components which are induced 
because of gravity, hydrodynamic and propulsion effects. 

K = Kp + K„ + K 1-13 

u n p 

Kg = (y^ cos 0 cos (j) - Zg sin (J) )w I-llj. 


The gravity term is: 





The hydrodynamic term is: 


Kx.= K,p + K.v + K qr + K rw + K vq + K pw 
^ p^ V qr^ rw vq ^ pw*^ 

+ K , ,v|w| + K , .v|u|+ K , ,v|v|+ K , ,p|u| 
vlwl ' ' vl u| ' ' v|v| ' ' P|up' 

+ K .r|ul+ K .pipl + K,. . 
r|u| ' I PIPI ' ' dist 


1-15 


The coefficients of this expression have "been determined hy model tests 
in a towing tank. 

Propulsion effects result from propeller action, effector action and 
movement of mercury within the vehicle. For the time heing, these effects 

will remain lumped under the tern K . 

P 

By combining Equations 1-12 through 1-15 the general equation describing 
the vehicle in roll may be obtained. In order to simplify the resulting 
equation, the relative magnitude of each term was determined by subjecting 
the system to the following set of standard conditions, and discarding the 
less important terms. The standard conditions are; 


u=v=w=1.0 ft/sec 

• • • . 5 

u = V = w = .1 ft/sec 

p = q = r = .01 Rad/sec 

• • • 

p = q = r = .01 Rad/sec 

(ft =z Q = \jj = ,10 Radians 

The roll equation which results is: 

^^xx" ^p^P “ sin(|>+Kp I-l6 


When the cross-coupling between u and p is considered, and if the assvunption 
of small roll angles (sin ^ is made, the basic equation describing 










vehicle roll is: 


(I_ - K‘)(j) = - MgZg(|) + K, 


XX p 


p|u 
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where K_ includes any moment about the x- axis applied to the vehicle, 

O 



Tank Arrangement and Coordinate System 


Tank 

Coordinates 

Pump To 

X 

y 

z 

T-1 

- 3.425 

0 

3.10 

T-2, T-3 

T-2 

- 3.425 

2.21 

- 2.24 

T-3, T-1 

T-3 

- 3.425 

-2.21 

- 2.24 

T-1, T-2 


Total mercury mass 82.5 slugs 
Vehicle constants of interest: 


M = 

4363 slugs 

Zg = 0.1335 ft. 

It 

X 

M 

56,200 slugs-Ft? 

K = - 5,300 slugs-Ft? 

p 

K* = 

V 

94T slugs-Ft. 

^v |u| ~ slugs (forward) 

Vlwl 

= 687 slugs 

K , = - 2,000 slugs-ft. 

p |ui 


Pumping rate (w) is related to commanded rate by: 

« 

^ dw • • 

• 083 + w = w 

at c 

Maximum pumping rate is 56.2#/Sec. 
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APPENDIX II 


PARAMETER CALCULATIONS 

The calculations illustrated in this section are based on the tank 
geometry shown in Figure two. The side tank area is assumed to be 0.6 
square feet, and the density of mercury is assiuned to be 26.3 slugs per 
cubic foot. 

The tank moment is calcvilated from Equation 2.6: 


”Ty = 


A(S)r sin 6 ds. 


Inspection of Figure one indicates ds = rd3. 
Therefore: 


Hpy = 2 


3^ 


A(S)r^ sin3 dB 


.<♦53 


.6 R^ sin B dB 


{y 

+ n 


2 A r^sin B dB 

J ^ 


'.<♦5 3 


For case I, r = 3.^35 ft. R = 3.06 


V 


= 2 {.562 + n A 3.40} ft*" 
P 


II-l 


For Case II, the similar equation is: 


= 2 {.562 + n A 3.68 } ft** 


II-2 


The solutions of Equations II-l and II-2 are tabulated in Table I. 

The value of the tank effective length is calculated from Equation 2.9: 


aTsT 


ds 
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N 

> 

•d 

" 

.00503 f 

.00637 

.00785 

.00955 

.0113 

,0133 1 

Case ] 

Case II 

I 

II 

I 

II 

I 

II 

I 

II 

I 

II * 

1 

1.21 

1116 

1.23 

1.17 

1.26 

1.18 

1.28 

1.19 

1.31 

1.20 

1.35 

1.22 

2 

1.29 

1.20 

I. 3 U 

1.22 

1.37 

I. 2 I* 

l.l*U 

1.26 

1.50 

1.29 

1.56 

1.32 

3 

1.38 

1.2lt 

l.lt5 

1.27 

1.52 

1.30 

1.60 

1.33 

1.69 

1.37 

1.78 

1.1*2 

It 

1.1*6 

1.28 

1.55 

1.32 

1.65 

1.36 

1.76 

1.1*0 

1.88 

1.1*5 

2.01 

1.52 

5 

1.55- 

1.32 

1.66 

1.37 

1.79 

1 . 1*2 

1.92 

1.1*5 

2.07 

1.53 

2 . 2 I* 

1.62 

6 

1.63 

1.36 

1.77 

I.I 12 

1.90 

1.1*8 

2.08 

1.52 

2.26 

1.62 

2.1*5 

1.72 

T 

1.'72 

1.U8 

1.88 

1.1*7 

2.05 

I. 5 U 

2 . 2 I* 

1.59 

2.1*5 

1.70 

2.67 ■ 

1.82 

8 

1.80 

l.ltU 

1.98 

1.52 

2.18 

1.62 

2 . 1*0 

1.66 

2.61* 

1.78 

2.90 

1.92 


Table I 


TANK MOMENT ABOUT PITCH AXIS 


A = 

P 

.00503 f 

t 

.00637 


,00785 

.00955 

.0113 

.0133 

Case I 

Case II 

I 

II 

I 

II 

I 

II 

I 

II 

I 

II 

.07 

.10 

.09 

.12 

.11 

.15 

.13 

.18 

.16 

.22 

.18 

.26 

.10 

.20 

.13 

.25 

.16 

.30 

.20 

.37 

.23 

.ui* 

.28 

.51 

.21 

.29 

.26 

.37 

.31 

.1*5 

.39 

.55 

.1*6 

,6k 

. 51 * 

.75 

.28 

.39 

.35 

.1*9 

.1*3 

.60 

.52 

.71 

.61 

.85 

.73 

1.0 

.35 

.1*8 

.1*1* 

.61 

. 51 * 

.73 

.65 

.90 

.70 

l.Ol* 

.90 

1 . 21 * 

.1*1 

.58 • 

.52 

.72 

.63 

.87 

.79 

1.08 

.90 

1.2l* 

1.07 

1.1*6 

.1*8 

.67 

.61 

.85 

.75 

l.Ol* 

.89 

I. 2 I* 

l.Ol* 

1.1*6 

1.23 

1.70 

.55 

.77 

.69 

.96 

.85 

1.15 

l.Ol* 

1.1*0 

1.19 

1.61 

1.1*0 

1.90 


Table II 


TANK FREQUENCY (SQUARED) g/L 
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This equation may be written for Case I as: 


L = 


453 


R d3 + 




453 


L = {1.39 + —^} ft. 


11-3 


For Case II, the similar equation becomes: 


L = {1.39 + } ft. 

P 


II-U 


The results of Equations II-3 and II-U appear in Table II, in the form 
of the parameter g/L. 

The tank curvature constant is defined in Equation 2.10 


M" 


= 1 


r cos Y <is 


For Case I y - 0 ds = rd3 


rj*53 


1^2 


M'' = 2 d 3 + 2 r^d 3 

•'O 'j*53 


= 


M 


2- „ 


3 U .8 ft? Case I 

16.0 ft^ Case II 


II-5 

II -6 


The tank moment of inertia is defined in Equation 2.11 


= prX(S) ds 


For Case I. 


= 2p { 


,4 5 3 


R' 


Ay d3- 




•45 3 


r^n A d3 } 
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f 






= {I 4 O 8 + 3070 N Ap } slugs ft2 


II-7 


From this equation the range of is U23 ^ £ 931 for the Case I geo¬ 

metry. 

Similarly for Case II: 



{U 08 + 9^.7 n Ap 


r^ds } 


From Figure tvo, 


r = 


1 . 3 ^ 
sin B 


j r= {lt08 + 680 n A } II -8 

• • J. P 

And similarly, the range of is ^12^ Case II geometry. 

Coefficient Evaluation 

The coef-Picient A is determined from Equation 2.25. 

5 


I' - K' + Jm 

X p T 


Using the values for I^ and given in Appendix I, and the results of 
Equation II-7 and II- 8 , it is revealed that A^ is approximately 1.62 x 10 
(slug-feet^) for both case I and case II geometry. 

A similar solution of Equation 2.27, using, the maximimi and minimum 


-s 


are: 


Case I 
Case II 


for both 

cases 5 

vas 

obtained. 

The results 

.316 < 

A 

2 

< 

.332 

Sec"^ 

II-9 

.316 < 

A 

2 

< 

.320 

Sec"^ 

II»10 


The evaluation of A and A are self-evident. 

3 <* 











Sample Damping Calculations 


Force = pA {1.5 + f(£/D)) V^/2 
P 

Case I il =8 feet Kinematic viscosity of mercury 
max 

— 6 2/2 

D = .13 feet Assumed to ^be 1.3 x 10 feet /sec 

n = 1.0 


V in 
ft/sec 

RxlO" 

®f(smth) f(r'gh) 

Smooth 

1.5 + 

Rough 

tUM 

Smooth Rough 
Lbs.F Lhs. 

Smooth 

F a 

Rough 

a 

1 . 

0.1 

0.0200 

o.o 44 o 

2.73 

4.60 

.478 

.735 

.478 

.735 

2 . 

0.2 

0.0175 

0.0320 

2.58 

3.47 

1.86 

2.51 

.465 

.627 

3 . 

0.3 

0.0160 

0.0290 

2.48 

3.29 

3.90 

5.17 

.433 

.575 

k. 

O.k 

0.0150 

0.0260 

2.42 

3.1 

6.77 

8.67 

.422 

.542 

5 . 

0.5 

O.OIU 5 

0.0250 

2.39 

3.04 

10.8 

13.3 

.432 

.532 

6 . 

0.6 

o.oiUo 

0.0245 

2.36 

3.00 

14.9 

18.9 

.413 

.524 

7 . 

0.7 

0.0135 

0.0240 

2.33 

2.97 

19.9 

25.4 

. 4 o 6 

.518 

8 . 

0.8 

0.0130 

0.0235 

2.30 

2.95 

25.8 

33.0 

.403 

.516 

9. 

0.9 

0.0128 

0.0230 

2.29 

2.91 

32.5 

41.3 

.401 

.506 

10 . 

1.0 

0.0127 

0.0225 

2.28 

2.88 

4 o.o 

50.5 

.40 

.505 

11 . 

1.1 

0.0125 

0.0220 

2.27 

2.85 

48.6 

61.2 

.402 

.506 

12 . 

1.2 

0 . 012 »t 

0.0215 

2.26 

2.82 

57.0 

71.0 

.397 

.492 

13 . 

1.3 

0.0123 

0.0210 

2.255 

2.79 

66.5 

85.0 

.394 

.502 

ll». 

l.U 

0.0122 

0.0210 

2.25 

2.79 

77.5 

95.6 

.395 

00 

CO 

• 

15. 

1.5 

0 . 0.22 

0.0210 

2.25 

2.79 

87.3 

108.0 

.388 

.480 

15. 

1.6 

0.0122 

0.0210 

2.25 

TABLE 

2.79 

III 

101.0 

124.5 

a 

ave 

.393 

. 4 l 4 

.486 

.533 


Sample Damping Term Calculations 

% " (smooth pipe) 

^ 1085 / V. • 

= -(rough pipe ) 
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The tank damping term is then defined; 


D = (smooth pipe) 

p 2pgA^ n ' ^ ^ ' 


- ~~ (rough pipe) 


Table III shovs the relation between D and A for N equal 1. 

P P 



A = ,00503 ft^ 

_ T) 

.00637 

. 00785 

.00955 

.0113 

.0133 

Case 

Rough 

Smooth 

Rough Smooth 

Rough Smooth 

Rough Smooth 

Rough Smooth 

R'gh 

Sm'th 

I 

3.92 

2.68 

2.71 

1.91 

1.85 1.37 

1.62 

1.20 

1.30 

.98 

1.07 

.83 

II 

3.01 

2.23 

2.23 

1.71 

1 . 7 U 1.33 

1.37 

1.07 

1.11 

.86 

.90 

.72 


Table IV 

Tank Damping As a Function of Pipe Area 


t 
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